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Abstract 

In this paper we obtain general integral formulas for probabilities in the 
asymmetric simple exclusion process (ASEP) on the integer lattice Z with near- 
est neighbor hopping rates p to the right and q = 1— p to the left. For the most 
part we consider an A-particle system but for certain of these formulas we can 
take the N — > oo limit. First we obtain, for the A-particle system, a formula for 
the probability of a configuration at time t, given the initial configuration. For 
this we use Bethe Ansatz ideas to solve the master equation, extending a result 
of Schiitz for the case N = 2. The main results of the paper, derived from 
this, are integral formulas for the probability, for given initial configuration, 
that the mth left-most particle is at x at time t. In one of these formulas we 
can take the N — ► oo limit, and it gives the probability for an infinite system 
where the initial configuration is bounded on one side. For the special case of 
the totally asymmetric simple exclusion process (TASEP) our formulas reduce 
to the known ones. 

I. Introduction 

Since its introduction nearly forty years ago (18] , the asymmetric simple exclusion 
process (ASEP) has become the "default stochastic model for transport phenomena" 
[21] . Recall [H [9] that the ASEP on the integer lattice Z is a continuous time Markov 
process r\ t where r]t(x) = 1 if x 6 Z is occupied at time t, and r) t (x) = if x is vacant 
at time t. Particles move on Z according to two rules: (1) A particle at x waits an 
exponential time with parameter one, and then chooses y with probability p(x,y); 
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(2) If y is vacant at that time it moves to y, while if y is occupied it remains at x. 
The adjective "simple" refers to the fact that the allowed jumps are only one step to 
the right, p(x, x + 1) = p, or one step to the left, p(x, x — 1) = q = 1 — p. The totally 
asymmetric simple exclusion process (TASEP) allows jumps only to the right, so that 
p = 1. 

In a major breakthrough, Johansson [6], building on earlier work of Baik, Deift, 
and Johansson [3J, related a probability in TASEP to a probability in random matrix 
theory. Specifically, if the initial configuration in TASEP is Z~ then the probability 
that a particle initially at — m moves at least n steps to the right in time t equals 
the probability distribution of the largest eigenvalue in a (unitary) Laguerre random 
matrix ensemble. The realization [6l [12] that the TASEP is a determinantal process 
[SI EE] has led to considerable progress in our understanding of the one-dimensional 
TASEP. (See [XT] for a recent review.) 

It is natural to ask to what extent these results for TASEP can be extended 
to ASEP. Since we no longer have the determinantal structure that is present in 
TASEP, random matrix theory methods, RSK-type bijections, or nonintersecting path 
techniques are not applicable (or at least not obviously so) to ASEP. 

However, it has been known for some time [U H] that the generator of ASEP is a 
similarity transformation of the quantum spin chain Hamiltonian known as the XXZ 
model [TH [20]. Since the XXZ Hamiltonian is diagonalizable by the Bethe Ansatz 
[20] . it is reasonable to expect that these ideas are useful for ASEP. Indeed, Gwa 
and Spohn [4] applied Bethe Ansatz methods to the TASEP for a finite number of 
particles with periodic boundary conditions (i.e., for particles on a circle) to compute 
the dynamical scaling exponent. 

Subsequently for TASEP on Z for a finite number of particles, Schiitz [15] showed 
that the probability that at time t the system is in configuration X = {xi, . . . ,x^}, 
given that its initial configuration was Y = {yi, . . . , y^}, is expressible as an N x N 
determinant. From this determinant representation Rakos and Schiitz [T3] derived 
Johansson's result relating TASEP to the Laguerre ensemble. The Rakos-Schiitz 
derivation uses the crucial fact that for TASEP the probability for any particle de- 
pends only on the initial positions for that particle and those to its right, and so it is 
expressible in terms of probabilities for finite systems when the initial configuration 
is Z _ . This is clearly no longer the case for ASEP. 

In this paper we obtain general integral formulas for probabilities in ASEP. For 
the most part we consider an iV-particle system but for certain of these formulas we 
can take the N —>■ oo limit, so that there are analogous formulas (involving infinite 
series) for infinite systems where the initial configuration is 

Y = {2/1,2/2, •• •}, 2/i < 2/2 < > +oo. 
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To specialize to infinite systems in TASEP with initial configuration Z we would 
replace Z~ by Y = Z + and let q = 1. 

Denote by Py(X; t) the probability that a system with initial configuration Y is in 
configuration X at time t. Then (Theorem 2.1) for an iV-particle system Py(X;t) is 
equal to a sum of N\ N- dimensional contour integrals. The integrand was suggested 
by the Bethe Ansatz (but there is no Ansatz!). The case N = 2 was established by 
Schiitz [15] in different notation, and he also proposed that there was a general result 
such as this, with some contours. 

The main objective of the paper is to obtain probabilities for the individual par- 
ticles at time t. To state the formulas we introduce some notation. We set 

e(0=pr 1 +ge-i, 

and define an A^-dimensional integrand I(x, Y,£), with variable £ = (£i, . . . , £jv) and 
parameter x G Z, by 

The parameter t appears in the last factor but not in the notation. 

For the probability for X\(t), the position of the first particle at time t, the sum 
of AH integrals given in Theorem 2.1 miraculously collapses into one integral. When 
we have 

nXl{t)=x)=p N(N-i)/2 f .../" I(XjY ,Od^---d^ N , (1.2) 

where C r is a circle centered at zero with radius r so small that all the poles of the 
integrand lie outside C r . 

In order to take the N oo limit we need integrals over large contours rather 
than small ones. This is because of the factors in the integrand fll.il) . Given 
a set S C {1, . . . , A^} we denote by I(x,Ys,£) the integrand analogous to I(x,Y,£) 
where only the variables & with i e S occur. Then when q ^ we have 

F(xi(t)=x) = y %2c s f ■ ■ [ I(x,Y s ,Od l % (1.3) 

where R is so large that all the poles of the integrand lie inside Cr. The sum runs 
over all nonempty subsets S of {1, ... , A^} and Cs are certain constants involving S 
and powers of p and q. 

Once we have this we are able to compute the expected value E,(xi(t)) and to take 
the — > oo limit in (11.31) . If Y = {yi, y 2 , . . .} with y 1 < y 2 < ■ • • — > +oo then on the 
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right side of (11.31) we simply take the sum over all finite sets S C Z + , the resulting 
series being convergent. 

These results are contained in Section III. In Section IV we derive the analogue 
of (11.21) for the second-left particle. 

The main results of the paper are in Section V where we obtain the analogues of 
(ll.2j) and (II .3p for the general particle. The analogue of (11.21) has the form 

¥(x m (t)=x)= °s>m f ■•■/ I(x,Y s ,0d^, (1.4) 

\S\>N-m JCr JCv 

where we sum over all subsets S of {1, ... , N} with cardinality at least N — m + 1, 
and cs, m is another explicitly given constant. The analogue of (11.31) is 

F(x m (t) = x) = Y, c s,m I •••/ I(x,Y s ,£)d% (1.5) 

\S\>rn ^ Cr ^ Cr 

where we sum over all subsets S of {1, . . . , N} with cardinality at least m, and cs, m 
is yet another explicitly given constant. The latter representation allows us to take 
the iV — > oo limit and so obtain probabilities for infinite systems. When the initial 
configuration is Z + the sum over sets in (11.51) may be simplified to a sum over integers. 

Of course ( 11.21) and ( 11. 3D and the second-particle formula are special cases of these 
but we give their proofs first, so that we can introduce the new ingredients gently. 

The deduction of these formulas from Theorem 2.1 requires two algebraic identities 
which were discovered by computer computation of special cases. The first is 



sgncT ( II ( p + ~ 

cteSjv \i<j 



t c2 t3 cN-1 
?o-(2) S«r(3) Sa(4) ' ' ' ^a(N) 



(1 - £<t(1)£<t(2)£<t(3) ■ ■ ■C<j(AT)) • ■ ■ (1 - ^ CT (Ar-l)^ CT (AT))(l ~ £<r(N)) 

- N(N-l)/2 ^li<j^ ~ ^ n 6] 

p a a-e, ) ' 1 j 

where the sum is over all permutations a in the symmetric group Sn- We also use an 
equivalent version of this identity, 
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X 



s s n ° ( n ( p + Q&M&u) - &«) 



(Ccr(l) _ 1)(£<t(1)£<t(2) - 1) • • • (£<t(1)£<t(2) * * * £a(N) 



i) 



N(N-l)/2 llt<j' 



- 1) 



;i-7) 



obtained from (11.61) by interchanging p and g and letting ^ 
For the second identity we introduce the notation 



iV-i+1- 



[iV] 



p 



A? 



q 



N 



p-q 



and then 



[N]\ = [N] [N 



[1], 



N 
m 



[N]\ 



[m}\ [N - m]\ 



'U 



where we set [0]! = 1. (Note that [*] is q m ( N m ) times a g-binomial coefficient with 
q equal to our p/q. Hence the notation.) The identity is 



e n 

\S\=m ies 
1 j£S c 



P + gj^j ~ & 

& - 6 



N — 1 

m 



N 



3=1 



1.9) 



for iV > m + 1. The sum runs over all subsets S of {1, ... , N} with cardinality m, 
and S c denotes the complement of S in {1, ... , N}. 

The proofs of these identities will be given in the last section. 

II. Solution of the Master Equation 

We denote by Y = {y 1 , . . . , y N } with y% < ■ ■ • < Vn the initial configuration of the 
process and write X = {x\, . . . , xn} € % N ■ When x\ < ■ ■ ■ < xn then X represents 
a possible configuration of the system at a later time t. We denote by Py(X;t) the 
probability that the system is in configuration X at time t, given that it was initially 
in configuration Y. 

Given X = {xi, . . . , xj^} e Z N we set 

Xi~ — . . . , Xi-i, Xi + 1, Xj+i, . . . , xn}, X i = {xi, . . . , Xi-i, Xi — 1, Xj+i, . . . , xn}- 

The master equation for a function u on 7* N x M + is 

d 



dt 



u(X; t) = J2 (p«(^5 *) + Q "W; t) - u{X- 1) 

i=l 



(2.1) 



and the boundary conditions are, for z = 1, . . . , N — 1, 
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U{X\, . . . , Xi, X{ + 1, . . . , Xjy] t) 



= pu(xi, ...,Xi,Xi,.. .,x N ;t) + qu(x t , . . . , + 1, x { + 1, . . . , x N ; t). (2.2) 
The initial condition is 

u(X; 0) = <5y(X) when X\ < • ■ • < Xjq. 

The basic fact is that if u(X; t) satisfies the master equation, the boundary conditions, 
and the initial condition, then Py(X; t) = u(X; t) when x± < ■ ■ ■ < x^u 

Recall that an inversion in a permutation a is an ordered pair {<j(i), <?(j)} in which 
i < j and a{i) > cr(j). We define 

c p + q£ a £p - 

Oa/3 ■ 



and then 
We also set 



p + q£a£p - £3 

A a = Y\_{Sa/3 '■ {o!,{3} is an inversion in a}. 



For now we shall assume p ^ 0, so the A a are analytic at zero in all the variables. 
Here and later all differentials d£ incorporate the factor (2ni)^ 1 . 

Theorem 2.1. We have when p ^ 

P Y (X;t) = [ ■ ■ f K \{€^ (l) ~ 1 ^ £ ^ t d^---d^ (2.3) 



ctGS 

where C r is a circle centered at zero with radius r so small that all the poles of the 
integrand lie outside C r . 

Remark. For TASEP with p = 1 we have 

S a/J = -(l-£a)/(l-&)- 

Fix a. In A a the factor 1 — £ a occurs for each inversion of the form {a, (3} and 
(1 — £ Q ) _1 for each inversion of the form {(3, a}. If a = a(i) then the number of the 



1 The idea in Bethe Ansatz (see, e.g., [71 [THl HO] ) , applied to one-dimensional TV-particle quantum 
mechanical problems, is to represent the wave function as a linear combination of free particle 
eigenstates and to incorporate the effect of the potential as a set of N — 1 boundary conditions. The 
remarkable feature of models amendable to Bethe Ansatz is that the boundary conditions for N > 3 
introduce no more new conditions, with the result that (|2.2[) involves only consecutive particles. The 
application of Bethe Ansatz to the evolution equation (master equation) describing ASEP begins 
with Gwa and Spohn [3] with subsequent developments by Schiitz [15] . 
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former minus the number of the latter equals a(i) — i. The number of inversions is 
the number of transpositions whose product is a. These give 

A. = sgn < 7jj(l-^ (i) ) ff «- i . 

Now the integrand in (I2.3P factors and we obtain the sum of sgn a times 

II / (i-^o)^^.r w ~ lce(e ' w) * d ^{<) = Il/ (l-ay^e^-'e^dt, 

i r i r 

and so the sum over a equals 

det ( I (1 - £) i_< C Xi ~ Vj ' % e £i0t <%\ . 

J Or 

This is the determinant representation of Py{X\t) obtained in [15] . 
To prove the theorem we shall show three things: 

(a) The right side of (12. 3ft satisfies the master equation for all X e Z N . 

(b) The right side of (12.31) satisfies the boundary conditions for all X e 7* N . 

(c) The right side of ( 12. 3D satisfies the initial condition when x\ < ■ • • < %n ■ 
Proof of (a). This is clear once the last factor in ( 12. 31) is written as the exponential 

O f Ei £ (&(i)K 

Proof of (b). We shall show that the boundary condition is satisfied pointwise by 
the integrand. Let TjU denote a with the entries a(i) and a{i + 1) interchanged. The 
boundary conditions will be satisfied provided thata 

for all a. Let us see why this relation holds. Let a = <r(i), (3 = a(i + 1), and suppose 
a > (3. Then {a,/3} is an inversion for a but not for Tjcr, so S a p is a factor in A CT 
but not in A^a, and all other factors are the same. Therefore, using S a /3 Spa = 1, we 
have 

A-Ticr = S/3 a A a = S^i+l).^*) At- 

The same identity holds immediately if f3 > a, since {/3, a} is an inversion for TjO" 
but not for a. Thus, (b) is established. 

2 The Bethe Ansatz proposes an integrand as in (|2.3[) with some coefficients A„. The fact that 
(|2.2j) only involves consecutive particles implies that satisfying it only requires relations between the 
coefficients A s and Ax ia - A straightforward computation shows that the relations are as stated. 
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Proof of (c). The initial condition is satisfied by the summand in (12.31) coming from 
the identity permutation id. So what we have to show is that 



when x\ < ■ ■ ■ < x^. This is the heart of the matter. We write 1(a) for the integral 
corresponding to a, and prove a series of lemmas. We think of a also as the ordered 
iV-tuple (a(l),...,a(N)). 

Lemma 2.1. Suppose that a appears in position a — 1 in a and that the entries 
preceeding a are a — 2 of the numbers less than a. Then 1(a) = 0. 

Proof. It follows from the assumption that there is a unique inversion of the form 
{a, (3} and none of the form {5, a}. Therefore the factor S a p appears in A a but no 
other S a s or Ss a - The variables £ a and £g appear in the integrand as 



We make the substitution 



ac 

power of £p that now appears is 



so that T) runs over a circle of radius r N , and we integrate with respect to £3. The 



The extra —1 in the exponent is due to the fact that d£ a = drj/£p. Since a can 
appear in no other S factor, the only one that can introduce a pole inside C r in the 
£0 integration is S a p, which becomes 



The apparent simple pole at ^ = coming from the third summand in the numerator 
does not occur because the exponent in (12. 4p is positive since yp < y a and x a > x a -\. 
The denominator is bounded away from zero when £p is inside C r , since |^| < r and 
the second summand has absolute value 0(r 2 ), so there is no pole inside C r and the 
integral is zero. 

Lemma 2.2. Suppose that in the permutations a and a' the entry a appears to the 
left of two adjacent smaller entries (3 and 7, and that the permutations differ only by 
an interchange of /3 and 7. Then 1(a) + I(a') = 0. 
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Proof. The pairs {a, (3} and {0,7} are inversions for both a and cr', so S a p S ai is a 
factor in both A a and A a >. Suppose for definiteness that f3 is to the left of 7 in cr, so 
a = cr(i), (3 = cr(j), 7 = cr(j + 1) with i < j. We make the substitution 



Then the powers of £g and £ 7 appear as 

gXj-Xi+ya-yg-l £Xj + i-Xi+y a -y^-l ^ 5) 

Both exponents are positive. We'll integrate with respect to £3 and £ 7 . Any S other 
than S a p or S a7 involving (3 or 7 will not introduce poles, as in the proof of Lemma 2.1. 
The product S 1 ^ becomes 

P + g *7 IL**,/? C 1 - C/3 P + g ?7 n^a, 7 ~ ^7 

We do the £g integration first. The first factor has a pole at zero as before but, as 
before, does not introduce one in the integrand because of the power of £g in (12. 5p . 
Assume for the moment that q ^ 0. Then the second factor has a pole at 

*>=rh n sr 1 - < 2 -6) 

The new second factor, its residue, has a pole of order 1 at £ 7 = (but again that 
does not contribute) and nowhere else. But the new first factor has a pole where 

n «.-'. p.7) 

with £g satisfying (12. 6p . Then £ 7 satisfies a quadratic equation with one of the roots 
inside C r , in fact 0(r 3 ), but all we will need is that (I2.6P and (12. 7p imply that £g = £ 7 . 

Now we compare this with cr'. The variables £g and £ 7 occur in different positions 
in cr and cr' and so their powers are different. After the variable change, for cr they 
are given by (12. 5p whereas for cr' they are 

eXj-Xi+ya— y-f — 1 £%i+i—Xi+ya— yp — l 

When we eventually have £g = £ 7 these are the same. The only difference, then, is 
the factor Spj. It occurs for a and not cr' when f3 > 7, and for cr' and not cr when 
/3 < 7. (Recall that /3 occurs to the left of 7 in cr.) It equals —1 when £g = £ 7 , so the 
sum of integrals equals zero. 
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Recall that we assumed q ^ 0. If q = then the integrals with respect to £p, after 
the substitution, for both a and a' are zero, so 1(a) = I(<j') = then. 

Lemma 2.3. For n > 1 the permutations in § ra can be grouped in pairs in such a way 
that the permutations in a pair differ only by one interchange of adjacent elements. 

Proof. Let two permutations form a pair if they differ by an interchange of the first 
two entries. 

Lemma 2.4. For any N the set §jv\{id} is the union of disjoint subsets, each of 
which consists of either of a single permutation satisfying Lemma 2.1 or a pair of 
permutations satisfying Lemma 2.2. 

Proof. We use induction, so we assume the result holds for N — 1. It clearly holds 
for N = 2, so we assume N > 2. For the set of permutations in which N appears in 
slot N we apply the induction hypothesis. Those permutations in which iV appears 
in slot N — 1 all satisfy Lemma 2.1 with a = N. Consider all those permutations that 
begin with a fixed (ai ■ • • a.N- n ~i N) with n > 2, so iV is in slot N — n. There 
are n\ of them, corresponding to the permutations of the n remaining entries. Pair 
these permutations as in Lemma 2.3. The pairs of permutations in Sn corresponding 
to these satisfy Lemma 2.2 with a = N. Putting all these together gives the desired 
decomposition of SV\{id}. 

Combining Lemmas 2.1, 2.2, and 2.4 completes the proof of (c), and so of Theo- 
rem 2.1. 

Remark. In case q ^ the same formula (12.31) holds when the circle is sufficiently 
large instead of small. A similar argument to the one just given should hold, but 
there is another way to see this. If we set 



and denote by P the probability density for the process with p and q interchanged, 
then P Y (X;t) = P Y -(X~;t). If we apply (12. 3p to this other process and then make 
the substitutions & — > £~ in the integrals we obtain (12. 3p with a large Cr. This 
duality will be use again in Section V, in the derivation of ( 11.51) from (11. 4ft . 



Here we determine the probability that the left-most particle x\ is at x at time t. 
Theorem 3.1. With C r as before and I(x, Y,£) given by (11.11) . we have when p ^ 



X 



{— xjv, ■ ■ ■ , — xi}, y 



{-y N , . . . , -yi) 



III. The Left-most Particle 



P(xi(t) =x)=p 



N(N-l)/2 



[ ■■■ [ i(x,y,o 



(3.1) 
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Proof. Since x\ < ■ ■ ■ < xn we may rewrite X = {xi, X2, ... , xn} as 

x, x + Zi, x + z\ + z 2 , ■ ■ ■ , x + z 1 H h Zjv-i- 

Then = x) equals the sum of Py{X;t) over all Zj, > 0. After summing, the 

integrand in f ]2. 31) becomes 

^ ( X ~ • • • ^Wj ^(2)C(3) • • • C(JV) TT Ux~ yi -l e e^)t\ 

" (1 ~ 6x(1)£<t(2) • ■ •C(T(Ar))(l - f CT ( 2 ) ■ ' £,a(N) ) • • • (1 ~ £ct(2V)) A A 

If we observe that 

YliP + Q£,<r(i)€*(j) - &(*)) 
A CT = sgn(T }| - w _ = sgncr — , (3.2) 



<r(i)>o-(j) 

then we see that the theorem follows from (jl.6p . 

We derive here the alternative expression for = x). Given a set 5 C 

{1,...,N} there is a corresponding set F$ = {yi '■ i & S} and corresponding \S\- 
dimensional integrand I(x,Ys,£). We define cr(S) = ^2 ieS i, the sum of the indices 
in S. 

Theorem 3.2. When q ^ we have 

a(S)-\S\ r r 

p N^^E j^» / ■•■/ J(*,n,0d |5| e, (3.3) 

where R is so large that all the poles of the integrand lie inside Cr. The sum runs 
over all nonempty subsets S of {1, ... , N}. 

We begin with a lemma that replaces integrals such as appear in (13.11) by sums of 
integrals over large contours. 

Suppose . . . , £at) is analytic for all & ^ and that for i > 

«-»•(£* -p)M* 

as — > 0, uniformly when all £j with j ^ k are bounded and bounded away from 
zero. Define 

t (t\ _ TT ^ ~ ^ /(^b 1111 6v) 

n,.(i-6) ■ 

For a subset 5 of {1,...,N} let If,s(£) denote the analogous function where the 
variables are the £j with i e 5, and in . . . , £/v) the £j with i £ S c are replaced 
by 1. 
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Lemma 3.1. Under the stated assumptions on . . . , we have when p, q ^ 

iri ^ = E .in^Z/s £ ■ ■ ■ £ rflsl « • < 3 - 4 > 

where r is so small that the poles of the integrand on the left lie outside C r and R is so 
large that the poles of the integrand on the right lie inside Cr. The sum runs over all 
subsets S of {1, ... , N}. When S is empty the integral is interpreted as /(l, . . . , 1). 

Proof. We use induction. The result is easily seen to be true when N = 1, so we 
assume N > 1 and that the lemma holds for N—l. We expand the £jv-contour on the 
left side. In addition to the pole at ^ = 1 we encounter poles at £jv = (£& — p)/q£k- 
We claim that the residue at this pole, when integrated over will give zero. The 
factor . . . , £jv), after substituting for £ N its value at the pole, is O(^) as — > 
by the assumption on /, while 

rL(i-&) ' 

after substituting for £jv hs value at the pole, will be of the order $,^ N+2 at = 0. 
The residue of l/{p+q£,k£,N—£,k) at the pole is l/q£k- The factor l/(p+q£i£,N— &) with 
i 7^ equals (p(Cfc — 6))- The factor — & is cancelled by the same factor in the 
numerator, so no new poles in the variable are introduced. So the l/(p + <z£i£jv — &) 
combined, including the contribution of the residue, give the power ef" 3 . Thus the 
product of all factors combined is 0(1). Hence the integral equals zero, as claimed. 

So after we expand the c^v-contour we have an integral where £jv is over an arbi- 
trarily large contour Cr and the other £j over small contours C r , and another integral 
(coming from the pole at — 1) m which ^ does not appear and the other £j are 
over C r . 

Let us consider the latter. The integral we get is the left side of (13.41) with N 
replaced by iV — 1 times 

R p +%zn-zL=i = ^- (3 - 5) 



i<N 



(Notice that £jy appears in the denominator of //(£) as 1 — C,n, and the pole at 
£n = 1 was outside the contour. The two minus signs cancel when we compute the 
contribution from the pole.) Our induction hypothesis tell us that this equals 

1 ^ 



V N-1 Q a(S)-\S\(\S\+l)/2 . . 

F SC{1,...,JV-1} H JCr JCr 

Now S c here indicates complement with respect to { 1 , . . . , N — 1 } . If we want to 
express this in terms of the complement in {1, . . . ,N} as in the statement of the 



12 



lemma we have to make the substitutions 

\S C \^\S C \-1, a{S c ) -> cx{S c ) - N, 
and therefore in the notation of the lemma the above equals 

Z2 Q a(S)-\S\(\S\ + l)/2 - " / J /^) dlS ^- 

SC{1,...,N-1} * JCr JLr 

This is the portion of the right side of ( 13. 4ft corresponding to those S not containing N. 

Now for the original integral, but where £/v is taken over Cr. If we integrate first 
with resepect to leaving us with an N — 1-dimensional integral with small C r , we 
see that we are in the case N — 1 with . . . , £n) replaced by 

m,...^..,)- [ n Jv/ 1 -p m :"f N) d ^. 0.6) 

We have to show that / satisfies the required condition. For notational convenience 
we take k = 1 and % = 2, so we make the substitution £ 2 — > — p)/q£i- All but the 
product are 0(£i) as £i — > 0, by assumption on /. Write the product as the product 
over % ^ 1,2, which is bounded if R was chosen large enough, uniformly for the £j 
bounded away from zero, times 

&v-fi Cat - (Ci ~P) hii 



p + <i€n€i - 6. p + Cjv (6 - p)/6 - (fi - p)Mi 



p + ?&v6 -6 p + £w (& - p)/6 - (fi - p)Mi 

The first factor equals while the second factor is Thus (13. 6p satisfies the 

required condition and we may again apply the induction hypothesis. 

If S C {1, . . . , N — 1} to compute § we replace the & with j € 5 C by 1 in 
/(Cij • • • 5 6v-i)- We see that the product in the integrand in (I3.6P is replaced by 



\S-\ ±± p+q £.£ N -£. 

ies 

If we set S = S U {iV} then in terms of S the full integrand including the £jv-variable 
is 



q\s 
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For the coefficient on the right side of (13. 4p with S replaced by S, the power of p is 
unchanged while the power of q is 

a(S) - \S\ (\S\ + l)/2 = a(S) -N- (\S\ - 1) \S\/2, 

and if we add to this \S C \ = N — \S\ we get 

a(S)-\S\(\S\ + l)/2, 

which is the power of q in (|3.4|) . 

Summing over these S gives the portion of the right side of (13 Ah corresponding 
to those S containing N, and this completes the proof of Lemma 3.1. 

Proof of Theorem 3.2. First assume p ^ 0. Observe that I(x, Y, £) is //(£) as 
defined above with 

/&,..., 60= (i-n&) n(^r*" v(6)t ). ( 3 - 7 ) 

and /(x, Fs, £) is If t s(£)- We show that / as defined this way satisfies the hypothesis of 
the lemma. The exponential summands £(£») and combine when £j = p) / 
to give 

+ Q — 7 V t + K* ~~ 2 = K* + 7 !> 



£fc - p ?6 Cfe £fc - P 

which is analytic at = 0. The powers of ^ and ^ combine as 



Cfc ~P 



fe 



as— l/i-l 



as 6c 0, which is 0(£&) since > y^. So the hypothesis on / is satisfied. Since 
/(l, . . . , 1) = the sum may be taken over nonempty subsets S. Because of the factor 
pN(N-i)/2 j n f |3 x H we must multiply the factor in (13. 4p by this, resulting in the factor 
in (Q- 

We can remove the requirement p ^ by taking the p — > limit since the power 
of p is nonnegative and no pole tends to infinity as p — > 0. 

An immediate conclusion from Theorem 3.2 is that F(xi(t) = x) tends exponen- 
tially to zero as x — > — oo, and therefore so does F(xi(t) < x). Thus F(x\{t) > x) 
tends exponentially to 1. 

We write 
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Q( X ) = P N ^ 2 1 ■■( er 1 - 1 • • • sr*- 1 n 



& - k 
P + qtitj ~ & 



where r is small. Clearly — > as x — > +00 and 

P(a;i(f) =x) =Q(x) + 

It follows that Q(x) = F(xi(t) > x), and this tends exponentially to 1 as 1 
Therefore 



-00. 



E( Xl (t))= lim Vy[Q(y)-Q(y + l)]= lim 

rr — ► — rvi * * t — ► — rv 



^Q(j/) + (x-l)Q(x) 



.J/=X 



lim 



J2Q(y) + ( X -l) 



y=x 



Now 



e w = ^~ i)/2 / • • 7 s ?/ 1 • • • v\ 1 n 



6 - 6 



X 



n*(i - &) 



(3.8) 



If we apply the procedure of the last section to this integral we start by moving 
the £jv-contour out. The resulting integral, with £at over a very large contour, is 
exponentially small as x — > —00, even though the other contours are small. As before 
the residues at the poles £jy — (£fc ~ P)l Q^k integrate out to zero. The contribution 
from the pole £jv — 1 gives the same expression we started with but with N replaced 
by N — 1. But there is now also a pole at £jv — ' " £jv-i whose contribution is 

-/ ••• / i/>n-i{€i, ■ ■ ■ ,€n-i;vn - yi, ■ ■ ■ ,vn - vn-i) d£i • ■ -d^N-u 



where 



X 



n 



p + qUZi ■ ■ ■ 6V)- 1 -6 1 - 6 • • • 6v n*(i - 6) 



J ptEi^W+eWi-^)- 1 )]* 
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For the right side of (13.81) when N — 1, which is what we are left with at the end, 
we expand the contour, encounter a double pole at £1 = 1, and find that it equals 



(p — q) t + ?/i — x + 1 + exponentially small term. 

Therefore 



N-l 



E(ar 1 (t)) = {p-q)t+y!-^2 / "/ ^i^ 1 ' ' • ' ' %'+i~2/i> • • • ' Vi+i-Vi) d ^ ' ' ' d ^r 

The integral f c ip\ (£ ; z) d£ has an explicit expression in terms of Bessel functions 
I n (2t). Indeed, it equals 

1 1 2-1 

2pte- 2t [4_a(2t) + 1,(2*)] + (2* - l)p {- e- 2t I (2t) -- + e~ 2t ^ 7,(2*)}. 

i=i 

The integrals of the other z) are not so simple. 

We now show how to obtain the probability P(xi(t) = x) for a system with 
infinitely many particles^ with initial configuration 

Y= {yi,y2, ...}, yi<V2< ► +oo (3.9) 

when q ^ 0. In fact, it is very easy once we have Theorem 3.2. We just modify the 
right side so that the sum runs over all finite subsets of Z + . Because Cr may be taken 
arbitrarily large the resulting series converges, as we shall now show. 

Consider the various factors in the integrand in (13.11) . where the £j are replaced 
by the £ n with n 6 S. The —1 parts of the exponents of the £ n may be removed if we 
replace d£ n by d£ n /£ n , which is 0(1). Suppose that \S\ = k as above. The analogue 
of the factor 

i-£i---6v 

n 4 (i-eo 

is at most (R k + 1)/(R - l) k = 0{2 k ) for R large. The product Ui<j^i ~ Q is at 
most {2R) k{k - 1 ^ 2 . The denominator Ui<j(P + Q^j ~ 6) is at least {qR 2 /2) k{ - k ~ l)/2 . 
The product of the is at most R kx ~^2y n . The exponential factor is 0(e aRk ) for 
some a depending on t. So the integral is 



O {{a/R) k ^' 2 R k *~Zvn e aRk ) = R°^-^ 



II n 



3 It follows from the fact that ASEP is a Feller process [8] that the limit as N — * oo equals the 
probability for the infinite system. We thank Thomas Liggett [10] for explaining this fact to us. 
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if we take R > a. (Since R is fixed the factor e aRk can be incorporated into the 
RO(k) term.) Since J2y n > J2(Vi + n ~ 1) = a (S) + ^ (Vi ~ 1) the above is at most 
RO(k)-°(S)_ And since a ( S j > ^ + this is at most R-^s)+o(a(S)^) _ 

The external factor in (13. 3p is 

<r(S)-fc fc(fc+l)/2- ff (S) < q -*(S)_ 

It follows that if we take R > 1/q 2 then the integral times the external factor is at 
most R~ a{ - S ^ 2 . 

Now consider all sets S with a(S) = k. Since the largest i G S is at most k, the 
number of such sets is at most 2 k . Hence the sum of the absolute values of the terms 
of the infinite series is at most a constant times 

oo 

Y^2 k R~ k /\ 

k=l 

which is finite when R > 4. Thus we have shown convergence for all t. 

IV. The Second-Left Particle 

In this section we compute the probability F(x2(t) = x). It is somewhat more 
complicated than that given for F(xi(t) = x) in Theorem 3.1 and the proof introduces 
some new elements. 

We use the notation (11.11) . and for 1 < k < N we set Yk = Y\yk- 

Theorem 4.1. With the contours C r as in Theorem 3.1 we have when p ^ 

N—l N—l r (■ 

P( X2 (t)= X ) = -g^ Z± p (N-l)(N-2)/2 I.I y- ^ 

N , x k-1 f f 

+ P ^ N -^ Yl ( - j / • • • / Yk > o dN ~^ 

k = \ \P / J C r J C r 

Proof. We rewrite X = {x±, X2, ■ ■ ■ , xn} as 

x — V, X, X + Zi,..., X + Z\ H V Z N - 2 - 

Then F(x2(t) = x) equals the sum of Py{X; t) over all v > and Zj > 0. If we sum 
first over z 2 , ■ ■ ■ , z^-2 the product Yii CaU)"^ m P-3p becomes, rL^"^ -1 times 

?^(3) W(4) ' " ' W(jV) 
(1 - £x(3)6t(4) ■ • "^(iV)) • ■ ■ (1 - £t(7V)6t(7V-1))(1 - £a(iV)) 
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We now move the £ CT (i)-contour out beyond the unit circle, and we do not encounter 
any poles. Here is the reason. From the first part of (13.21) we see that we get poles at 

{(&(*) -p)/q&(k) if Ac > i, 
(4.2) 
p/(l-?6r(k)) iffc<i. 

The poles in the £ CT (i)- variable are at = (£ CT (fc) —p)/q£a(k), and these are very large 
since ^(ja G C r . So we move the ^(^-contour out beyond the unit circle and sum, 
giving 



(4.3) 



(4.4) 



1 C f 2 cN-2 

1 S<r(3) W(4) " " •S a (Af) 

£<r(l) - 1 (1 - £7(3)67(4) • "€<t(N)) • • • (1 - €tr(N)€(r(N-l)){l - 6<x(iV)) 

If we move the contour back to C r we pass a pole at £ a m = 1 with residue 

t tl cN-2 
£<x(3) ? CT (4) ' ' 'W(jV) 

(1 - €7(3)67(4) • ■ ' ^a(N) ) • • " (1 - £<T(A0£o-(iV-l))(l - £<r(JV)) 

The factor A a when = 1 equals 

sgna — IN , ' : 1 < h cr(i) > cr(j , i,j > H, (4.5) 

since there are <r(l) — 1 indices less than er(l). 

Now we add all terms in which cr(l) = k. For the contribution from the pole at 
Co-(i) = I; when <r(l) = (14.51) may be written (as in (13.21) ) 

sgna TT< 77 — 77 II(p + - 

where all indices are 7^ fc. This is to be multiplied by (14.41) . If we use identity (11.61) 
to sum the product over those a with cr(l) = k we get 



/ x 11(^-6 



4 / n^-^n^+^-ft 



J? 



where again all indices are 7^ k. The exterior factor is now Yltyk (£f V ' 1 e6 ^ t ), an d 
from these we obtain the sum on the right side of (14.11) . 
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Next we consider (j4.3j) . which we rewrite as 

(1 — £tr(2)£fr(3) ' ' ' £cr(N)) ^ <J ( 3 ) ^(4) ' ' ' £<j(jV) 



6c - 1 (1 - C<t(2)Co-(3) • ' £c(iV) ) • • • (1 - ^(AT)Co-(iV-l))(l - £<r(iV))' 

where a is now a map from {2, . . . , iV} to {1, ... , iV}\{A;}. The factor sgncr becomes 
(— l) fc+1 sgncr. We also rewrite A a as 



If we use identity (JL6J) to sum this over these a we get 

The other factors are still Y\ i e^» e ^*). 

To evaluate the sum of 04.61) over k we write it as 



-P 



(JV-l)(JV-2)/2 i^i 



no - n(p 



X 



The factor (— l) fc became —1 because of the way we rewrote the product of the £j — 
in (14.61) . Identity (11.91) with A; = 1 tells us that the last sum equals 



pN~l _ gN~l 

Q 

p-q 



i-n« 



J 

3 



If we recall the power of p in the first factor above and the ubiquitous factor Y\ { Vl 1 
we see that this gives the first term in (14.11) . 
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V. The General Particle 



In this section we consider the mth particle from the left for general m. We prove, 
with the notation (jl.8p . 

Theorem 5.1. We have when p^O 



x (-l) m -Hs c l 

\S c \<m 



Xm (t) = x) = p(N-m){N-m+l)/2 q m(m-l)/2 

r I(x,Y s ^)d l % 



\s\- 

m — \S C 



qa(S c )-m\S c 



pj(s°)-\s*\<\s°\+i)/2 



(5.1) 

where r is so small that the poles of the integrand lie outside C r . The sum is taken 
over all subsets S of {1, . . . , N} with \S C \ < m. 

We shall first establish a preliminary form of the result, and for that we use a 
lemma analogous to Lemma 3.1 and which follows from it. Now we have a function 
9 (£,!■> • • • ■> 6v) which is analytic for all & ^ and satisfies, for i < k, 

0(6..... fr) I /n , =0(C X ) 

as 6c — * 00 ? uniformly when all £j with j ^ k are bounded and bounded away from 
zero. Define, as before, 

j (t\ _ TT ~ ^ • • • ,6v) 

9lU lip 11,0-6) ' 

and similarly I Sj s(£). 

Lemma 5.1. Under the stated assumptions on ■ ■ ■ > 6v), we have when p, q 7^ 

,a(5 c )-7V|S c | 



/ - / j.(t)^{=D- 1 ) w pW g,4x tW )^ w / - 1 we^f. (5.2) 



where r is so small that the poles of the integrand on the right lie outside C r and R 
is so large that the poles of the integrand on the left lie inside Cr. As before, S runs 
over all subsets of {1, . . . , iV}. 

Proof. We apply Lemma 3.1 with p and q interchanged to the function 

/(a,---,6v)=#(G\...,er 1 )nc 1 - 

The required hypothesis on / follows from the hypothesis on g. The left side of ( 15. 2p . 
after the substitutions 6 — > l/6v-i+i ; equals the left side of (13 ,4p after interchanging 
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p and q, times (— 1)^, because Y[{^ ~ 6) becomes rite ~ !)• So we apply Lemma 3.1 
to / and then change variables again resulting in a factor (— l)' 5 ' in each summand. 
The reason for the different coefficients is the reversal of the order of the variables. 
Thus the coefficient of the integral involving If^s with p and q interchanged equals 
the coefficient of the integral involving I §, where S = {N — i + 1 : % G S}. This, 
together with the resulting power of —1, is what we see on the right side of (15. 2\i . 
This proves Lemma 5.1. 

To state the preliminary form of the result we introduce more notation. For 
disjoint subsets T and U of {1, ... , N}, we define 

I(x,Y Tt u,Z) 



n n (,, ■ qu, c: 



i-m ^Trrn ^nT+^ n lifer- 1 (5.3) 

i i<j 

where indices with unspecified range run over T U U. If T C U we define cr(T, U) to 
be the sum of the positions of the elements of T in U. Thus, if U = {2, 3, 5} and 
T = {2, 5} then a{T, U) = 1 + 3 = 4. Finally, for a set U we define 

sgaU = (-l)#{( i J) :i> J' i6 ' 7 'J eC/c }. 
Lemma 5.2. With r small enough and p, q ^ we have 

M Xm (t) = X) = p( N ~ m )( N ~ m + 1 )/' 2 + m (m-l)/2 ? (m-l)(m-2)/2 

_Xf/\T)-(m-l) |C/\T| 



\U\=m-l TcU 



p a(U\T)+\T\(m+\U\T\)/2 



X 



j(x,y TiC/c ,Orf |TUC/ \, 



where U runs over all subsets of {1, . . . , iV} with \U\ = m — 1. 
Proof. To begin we now write X as 

X-V m -i Vi, X-V m ^ 2 Vi, ••-,X-Vi, X, X + Zi,..., X + Z\-\ h^Ar- m . 

The product [7* fo)^" 1 in ^ is re P laced b ^ IL (^~ w-1 e e te)*) times 

— HI »m-l t~ Ul £«1 f^lH hZjV-m 

W(l) ' ' ' W(m-l) W(m+1) ' ' ' W(jV) ' 
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and we have to sum over all V{ > 0, Zi > 0. As in the proof of Theorem 4.1 we can 
sum over the immediately and we get 



-vi Vm-l £-vx W™+l)W(m+2) ' ' ' W(jy) 



J-—V1 V m —1 J-- 

W(T) ' ' ' ?(T 



(1) CT(m ^ (1 - £<r(m+l) ■ • ■ ia(N)) " • • (1 - ia{N)) 

Before we can sum over the Vi we have to move the £ CT (j)-contours out, and we do 
them in the order i = 1, . . . , m — 1. As in the proof of Theorem 4.1 we see by referring 
to (14.21) that the poles obtained from moving the £ CT (i)-contour are very large and so we 
can move that contour out almost that far. Then if we want to move the ^(2) -contour 
out we encounter poles with k > 2 in (14.21) . which are far out and so no problem, but 
also the pole with k = 1 when <x(2) < ct(1), which is at £ CT ( 2 ) — p/(l — q(,a(i))- We 
show that the residue at this pole, when integrated with respect to £o-(i), gives zero. 

Recall that the ^ CT (i)-contour is large, and so £ CT ( 2 ) as a function of is analytic 
outside the £ CT (i)-contour and is in fact 0(£~^) at infinity. The part of the residue 
that comes from the factor 

P + g&r(l)£<r(2) — £<r(l) 
P + gCa(l)C<x(2) - £<t(2) 

is 0(1) at infinity as are all the other factors in the first part of (13.21) because £ CT (2), 
in terms of , is small when ^ ff ^ is large. The powers of the & involving and 
£ CT (2) combine as 

/ \ -vi v m -2-x-l f r> \ ~~ y " { V 

which is analytic outside the ^ cr ( 1 )-contour and O^^' 1 VaW 2/ct(2) ) at infinity. The 
exponent is < —2 since > and y a {2) < Va(i)- Finally we have to check the 

exponential of J2 £ (£i)- The sum of those involving ^(i) and £ CT (2) is 

y— + + (1 - <?£x(i)) + q- ^7 , 

which is bounded at infinity. Hence the -integral is zero, as claimed. 

Continuing this way we move all the £ CT (j)-contours out for i < m. We then sum 
over the Vi, obtaining 

1 



(6<t(1) - 1) (6t(1)£<t(2) - 1) • • • (£ CT (i)£ CT (2) • • -C<T(m-l) ~ 1) 

e c2 £N—m 

?o-(m+l)? -( m+ 2) ' ' 'S cr (AT) 



X- 



(1 - £<r(m+l) • ■ ■ £<r(AT)) ■ ■ ■ (1 - £a(N)) 
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as replacement for what we had before. We also have, from the numerator in the 
second part of (13.21) . a product which we write as 

n (p + s&w&g) - &«) n (p + q&w&u) - n (p + - 

i<j<m m<i<j i<m<j 

Therefore we are to take the sum over all cr of sgn a times the product 

1 



II (P + - &(<)) 77 7777 7" 



,j(2) - 1) ••• (6r(l)&(2) ••• £<r(m-l) - 1) 

t t2 tN—m 

?cr(m+2) i ; tT ( m+ 3) ' ' ■So-(JV) 



[P + QKa{i)Ka(j) - 4<x(i)) — -, t \ n — 7 — 



X 

m<i<j 

X IJ (P + q€*(i)£a(j) ~ 6r(t) 
i<m<j 



Now we take a fixed [/ C {1, . . . , iV} with |/7| = m — 1, and sum over all permu- 
tations a such that a(i) £ U when i < m — 1. Notice that the last product above is 
equal to 

H(p + qtej-ti, (5-4) 



i6£/ 



and so is independent of the particular permutation. Therefore we may sum inde- 
pendently over bijective maps {1, . . . , m — 1} — > U for the first factors and bijective 
maps {m, . . . , N} — > U c for the second factors. 

To keep track of the signs of the permutations we observe that sgn a equals the 
product of the signs of the two restrictions of a times sgn U. So we eventually have 
to multiply by sgn U. 

For the second factor above we use (11.61) and find that the sum equals 



n - 6 



p (N-m)(N-m+l)/2 



I'' ) II 



i€U c 

i€U c 



For the first factor we use (11.61) and find that the sum equals 

n & - 6) 



(m-l)(m-2)/2 (5 Q) 

rife- 1 )' 
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This is to be multiplied by (15. 4p . divided by rii<j(P+s£i£i — &)> an d then integrated 
over large contours Cr with respect to the & with i E U . But we want all integrals 
to be taken over the same contours C r so we want to replace the integral with all 
contours Cr with a sum of integrals with all contours C r . 

In our application of Lemma 5.1 N will be replaced by m — 1, {1, ... , N} will 
replaced by U, and S will be replaced by T, which is the reason we defined cr(T, U) 
as the sum of the positions of the elements of T in U. The coefficients become in this 
notation 



q a(U\T, U)-(m-l)\U\T\ 

p \T\(m+\U\T\)/2-a(T,U)- ( 5 ' 7 ) 

If we put all integrands together the result is, aside from the powers of p and q in 
5SD and (1531) . 



i<3 



i-n&)n&-&) 

n e e ^') (all indices in f/ c ) (5i 



i 

11(0 - 6) 



X 



g + g&O - 6 
p + - & 



n (^r yi_1 e£te) ') ( a11 indices in f/ ) ( 5 - 9 ) 



* n 

We apply Lemma 5.1 with 



i>3 

ieu,jeu< 



i 6 (7 •'U- "'w i>:j f ^s«sj - 6f/c 

i£U,jeU c 

where gr(£) = 0({&}i 6 ff). 

The poles of the integrand defining g are at 

& = K, - /'). <,• 

where j G t/ c . Since £j can be arbitrarily small the pole is outside Cr for any R, so g 
is analytic for all £j 7^ 0. 
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To check the main hypothesis on g we observe that after the substitution £j — > 
pj (1 — qik) the product C[~ Vi ~ l £fc~ yfe_1 becomes 0{^{!~ yk ) as — > oo which is 0(^ A T 1 ) 
since i < k. The sum of the exponents in the last factor which involve and is 

it Pit 

which is bounded at infinity. As for the integral, we show that it is bounded when 
one £fc — > 00 while the others are bounded. If we take a fixed j > k with j G U c , the 
pole at £j — p/(l — q^k) passes across the ^-contour when ^ — > 00. The residue of 
the factor in the product that contributes the pole is seen to be 0(1) as are the other 
factors involving % 7^ k. We do this with each j in turn and end up with a sum of 
integrals each of which is 0(1). If we make the substitution & — > p/(l — q£k) then 
£fc — » 00 while — q^) remains bounded, so g satisfies the required hypothesis. 

To find the summand in (I5.2p we must evaluate g(£) where all the with i G U\T 
set equal to 1. Each factor in the product in the integrand with such an i is (—q/p). 
For each such i the number of j in the product in the integrand satisfying i > j and 
j G U c equals % minus the position of % in U. The sum of this over all % G U\T equals 
a(U\T) — a(U\T, U). If we multiply (I5.7P by (—q/p) to this power the result may 
be written 

n a(U\T)-(m-l) \U\T\ 
(_-q<r(E7\T)-«r(tf\T,lO ( / 



pcr(U\T)+\T\(m+\U\T\)/2-m (m-l)/2 ' 



since a(U\T, U) + a(T, U) = m{m - l)/2. 

There are also the factors (— i)\ u \ T \ coming from (j5.2p and (— l)' 17 ' coming from 
the fact that Yl(d ~ 1) appears in (15.61) rather than FJ(1 — ^). These factor combine 
as (-1)I T L 

For the integrand we must combine (15.81) . (15.91) with the ^ with i G U\T set equal 
to 1, and the integrand in g with the £j with z G Z7\T set equal to 1. The result is 
(15. 3p with U replaced by U c . If we take account of the original factors in ( 15. 6ft and 
( 15. 5p we have established Lemma 5.2. 

Proof of Theorem 5.1. Suppose temporarily that q 7^ also. We take a fixed 
S* C {1, . . . , iV} with \S C \ < m and in Lemma 5.2 sum over all T and U with T C U 
and T U U c = S. Let us write everything in terms of T and U c . 

For any i E U, the position of i in [/ equals #{j : j < i, j G U}, so 
<r(T,U)=#{{i,j):i>j, ieT,jEU}. 

In particular, 

a (U\T, U) = #{(z,j) :i>j,ie U\T, j G [/}, 
(7(C/\T) = : i > j, i G f/\T, j G {1, • ■ ■ , iV}}, 
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and so 

a{U\T) - a(U\T, U) = :i>j,ie U\T, j e U c }. 

Also, 
Thus 

(_l)|r|+#{(i,j):i>j,iGT,jG(7 c }_ (5.10) 

is the combined power of —1 that occurs. 

The powers of p and q that occur in the summation are, since U\T = S c , 



„cr(S c )-(m-l) \S C \ 
p a(S-)+\T\(m+\S-\)/2^ 



(5.11) 



and only depends on |T|, given S. 
In (15.31). we write 



n 

n & - = 



i,jeu c ur 



II n 



i,jSU c or i,j€T 

ieu c , jer ieT, jeu c 



The denominator may be written 

n n (o-6)=(-i) #{( ^ ieT - wc} n 

»<J iGT, 7'GC/ c 

The power of —1 combined with (I5.10p equals (— 1)' T ' and therefore our integrand, 
aside from this power of —1, may be written 

n n&-&) 

(i-TT f) ieT ' J&u ° i^i TT ft*-**- 1 P «(eo A 

1 ,L A ; n nc 1 - 6) n(p + ^ - ¥ ^ ;> 

i£T,j£U c i i<j 

where indices not specified range over T U U c . 

Now we take a fixed S C {1, . . . , N} with \S C \ < m and in Lemma 5.2 first sum 
over all T and U such that T U U c = S. The condition \U\ = m — 1 translates to 
|T| = m — 1 — |>S C |. Apply ( 11.9jl with {1, . . . , iV} replaced by S, with S replaced T, 
and with m replaced by m — 1. We obtain (15. ip . and this completes the proof when 
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We can remove this condition by taking the q — > limit since no pole tends to 
zero when q — > 0. 

We now obtain the expansion where all integrals are taken over large contours. 
Theorem 5.2. We have when q ^ 

¥{x m {t) = x) = (-1)^+1 ( pg )m(m-l)/2 

^-|5|(|5| + l)/2 / / '(^,0^, (5-12) 

JCr J Cr 

where R is so large that the poles of the integrand lie inside Cr. The sum is taken 
over all subsets S of {1, ... , N} with \S\ >m. 

Proof. Denote by P the probabilities for the process with p and q interchanged. 
As in the remark following the proof of Theorem 2.1, F(x m (t) = x) is equal to 
F(xN- m +i(t) = —x) with initial configuration Y replaced by {— y^, . . . , —yi}. In the 
integrals in (15.121) we make the replacements £j — > l/£/v-i+i- The upshot is that in 
( 15. ip we replace m by N — m+ 1, in the coefficients we replace S by S 1 , we multipy by 
(— 1)' 5 ' +1 (because of the sign change in the integrand), and take the integrals over 
Cr. A little algebra, using the general fact 

a(S) = -i + l) = (N + l)\S\- a(S), 

ies 

shows that the result is (I5.12p . This complete the proof of Theorem 5.2. 

As with the first particle when Y is infinite and bounded below, we can show that 
the sum (15.121) converges when it is taken over all finite subsets of Z + . This gives the 
probability for infinitely many particles. 

In the special case Y = Z + we can evaluate the sum over all sets of a given 
cardinality. We define the fc-dimensional integrand 

ij p + - ^ n( i - 6) - p) i 

i 

where all indices run over {1, . . . , k}. The result we obtain is 
Corollary. When Y = Z + we have when q ^ 

F(x m (t) =X) = (-l) m+1 q m(-m-l)/2 
p (k-m)(k-m+l)/2 q k(k+l)/2 I ... I J fc ( Xj ^ d ^ . .. d ^ (5.13) 

J Cr JCr 



|S|-1 
IS"! — m 



k 

k- 



- 1 

m 
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Proof. We sum first over all S C Z + with \S\ — k. If S = {z±, . . . , z^} we make the 
variable changes £ Z1 — > £1, . . . , £ Xh — > in the integral of J(x, Y5, £), so all integration 
are over the variables & with « e {1, . . . , The integrand becomes 

where all indices run over {1, . . . , fc}. The only part of the coefficient in f 1 5 . X 2 j) in- 
volving more than \S\ is (p/q) a ^ s \ When we multiply by this the product may be 
written 

rife - 6) 



X 



If we sum over all {zi} with < z\ < ■ ■ ■ < z^ the last product becomes 

l) ((j 6) • • ■ (; &) - 1) ((; 6) • ■ ■ (| &) - 1) • • • ((* &) - 1 

The integral is unchanged if we antisymmetrize this, because all other factors are 
symmetric except for the Vandermonde. We can bring this antisymmetrization to the 
form of identity ( 11.71) if we make the substitutions 

, P 

q,i — ~ TJk-i+i- 
Q 

The second factor becomes 

1 



{r]i - 1) (77! T]2 - 1) • • • (771 772 • • • rj k - 1) ' 
while the first factor (after the index changes % — > k — % + 1, j — > k — j + 1) becomes 

fc(fc-l)/2 

i<j 

Now we can apply (11.71) with p and q interchanged and we see that the antisym- 
metrization of this is 

_|_ _P 2l£ = _ D *(fc+i)/ 2 Hi 

kl qk{k -m Yl (vi _ 1) A;! P n^-P)' 
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If we recall the remaining factor 



p -m\S\ gk(k+l)/2 



g -|5|(|5|+l)/2 pmk 

in the coefficient in (15.121) we see that we obtain formula (I5.13p . 

Remark. The power of p on the right side of (15.131) is always nonnegative and is zero 
only when k = m. Hence when p = 0, in other words in the TASEP where particles 
move to the left, only one term survives and we obtain 

f_-iWm-l)/2 r r t . . . t _ i 



ml 



i 

For F(x m (t) < x) we sum ¥(x m (t) = y) over all y < x (which we may, since R > 1), 
obtaining 

(_1 )m(m-l)/2 r r 

n*m(t) < x) = — / •■•/ nfe - life- ^ 
xll(r ra ^ rl)t ) di 1 ---du. 

i 

By a general identity [2] this equals 

(_-Qm(m-l)/2 det ^ ^i+i+x-m ^ _ ^-m e (£-l)t 



i, j=0,...,m— 1 



After reversing the order of the columns this becomes a Toeplitz determinant equal 
to the determinant of Rakos-Schiitz [T3| (12)] which they used to obtain Johansson's 
result [6|. 



VI. Proofs of the Identities 



Proof of identity ( 11.6j) n We use induction on N, and assume that the identity 
holds for N — 1. (It clearly holds for N = 1.) Call the left side <p N (£i, ■ ■ .,£n). We 
first sum over all permutations such that er(l) = k, and then sum over k. If we 

4 When we showed Doron Zeilberger the identity when it was still a conjecture he suggested [22] 
that problem VII. 47 of [11], an identity of I. Schur, had a similar look about it and might be proved 
in a similar way. He was right. 



29 



observe that the inequality i < j becomes j ^ i when i — 1, we see that what we get 
for the left side is 

1 N 

— — ^(-l)^ 1 TUp + qZktj ~ 6) TT & ■ ^v-i(6, • • • , 6-i, 6+i, • • • , 6v), 

which may also be written 

: ^TT-^h EC" 1 )** 1 11^ + ~ &) ■ C 1 <P*-i(£i> • • • , 6-i, 6+i, • • • , 6v). 

We want to show that this equals the right side of (11.61) . and the induction hypothesis 
gives 

n 

^-i(6, • • • , 6-i, 6+i, ...&) = Z^ 1 ^/ 2 . 

ri( i - ^ 

After some multiplying, dividing, and computing powers of —1 we see that what we 
want to show is 

til 74 6 1 Ij&KSj ~ 6) 4i 42 • • • 4tv 

If we change the first product on the left to run over all j we have to divide by 
P + q£,k ~ 6- Setting p = 1 — q shows that 

1-6 i 



p + q£l ~ 6 p- ?6 ' 

So the left side of (16.11) equals 

N N x 1 

We evaluate this by integrating 

W 1 1 
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over a large circle. Since the integrand is 0(z 2 ) for large z the integral is zero. There 
are poles at and the and the sum of the residues there is 



N N 



There is also a pole at z = p/q and for the residue there we compute 



V + Vjj -p/q _ q + g£j - 1 

ij-p/q q^j-p 

so the residue at p/q is —p N_1 . This gives 

N N . 



P, 



P 



AT-1 



E 11/' + <** - 60 ■ ^(^r^y n .^ K ._ &) n . (i 



P 



JV-1 



as desired. 

This completes the proof of identity (11.61) . 
Proof of identity ( II. 9|) . We use the easily established recursion formula 

"N- ' 



N 
m 



P 



N-l 
m 



+ q 



N—m 



and a preliminary simpler identity, 



E n 

\S\=m *eS 



p + g&£ 3 - - 6 
- 6 



m 

N 
m 



(6.2) 



(6.3) 



where S is as before. 

We prove thi£@ by induction on N, so assume (16.31) holds for N — 1. The left side 
is symmetric in the & and is 0(1) as any & — ► oo with the other fixed. If we 
multiply the left side by the Vandermonde rL<j(£« — £j) we obtain an antisymmetric 
polynomial which is 0(^ J Ar ~ 1 ) as any £j — > oo with the others fixed, so it has degree 
at most iV — 1 in each £j. Being antisymmetric it is divisible by the Vandermonde, 
and having degree at most iV — 1 in each of the £j separately it must be a constant 
times the Vandermonde. Thus the left side of (16.31) is a constant, say Cjv, m - 

To evaluate the constant set £jv = 1. For convenience we write 



P + q££j ~ & 
" 6 



3 The proof is a modification of one found by Anne Schilling [TJ] for an equivalent identity. 
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We have 



c N , m = j2 n u &>ti)i _ = <i N ~ m e n ufati) +p m e n ^.^) 



5 «es 



NeS ies\iv 



jes c \jv 



By the induction hypothesis the right side equals 



N-m 



' N — 


1" 




' N — 1" 






+ p m 




m — 


1 


m 



and this equals [ 1 by (16.21) . This establishes identity (16.31) . 

The proof of (11. 9p runs along the same lines. We interpret both sides to be zero 
when N = m, and do an induction on N > m. So we assume N > m and that the 
formula holds for N — 1 > m — 1. We quickly deduce that the left side is a polynomial 
of degree at most one in each £j. If we call the left side CAr, m (£) then 



CW,m(Cl) • • • >£jV-l, 1) — <? 



N—m 



NeS ies\N 
jes c 



ih 



+P m E n n 



j£S c \N 



q 



j&S c \N 

Jto- \JV 

m Cat — l,m— 1 



7V-1; 



Similar relations hold for the other £j. Notice that when N = m the second sum above 
does not appear and C m _i )m _i(£) = so this is consistent with our initial condition. 

If we call the right side of (11.91) C' Nm (£) we see from (16.21) that the same re- 
lations hold for C' Nm (£) (with initial condition C' mm (^) = 0) and so for the dif- 
ference Z)jv,m(0 = CV,m(0 — C' Nm (^). The induction hypothesis gives D N _ lm = 
DN-i,m-i = 0, so we have shown that for any i 



D N ,m(0\ 



6=1 



0. 



Any polynomial which has degree at most one in each £j and vanishes when any £j = 1 
is a constant times Yli^i ~ so Aw,m(0 nas this form. 



We have shown that 



«)=«)+cnte-i) 



6 Such a polynomial must be of the form £jv — 1 times a polynomial in £i, 
property, so the statement follows by induction. 



, , £at_i with the same 
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for some c. We show that c = by computing asymptotics as £jv — > oo. All terms 
are asymptotically a constant times £at. If in the sum in (II .9p N £ S then the 
corresponding summand is 0(1). So we need consider only those S for which N (jL S. 
In the product in the summand, if j = N then the corresponding product over i has 
the limit q m Ylies& since there are m factors with limit q^. It follows that 

e ii E± #% zi -n«-- n «. 

S|=m iSS s -? s * i e s j<N,j£S c 

SC{1,...,JV-1} iSS c ,j<N 



,. CV m (£) 
hm — = — 



p + g&6 - ^ 



II E II 

i<JV |S|=m i6S Sl 

Sc{l,...,N-l} j£S c ,j<N 



Identity (16. 3p tells us that this equals 



n a 

i<N 



N — 1 
m 



Clearly C^r m (^) has the same asymptotics, so c = 0. 
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